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Intermediate-frequency-peaked BL Lacertae objects (IBLs) are a class of blazars characterized 
by a spectral energy distribution (SED) with a lower-energy synchrotron peak than the majority 
of extragalactic sources detected by ground-based imaging atmospheric Cherenkov telescopes 
(lACTs). Consequently, the peak gamma-ray flux falls outside the very-high-energy regime 
(VHE, >100 GeV) covered by lACTs such as VERITAS, making IBLs difficult to detect except 
during infrequent episodes of elevated flux. However, the study of these sources in a low-flux 
state is essential for developing a complete understanding of the blazar paradigm. We present the 
results of VHE analyses of long-term low-flux-state observations completed for two IBL sources: 
3C 66A and W Comae. Eor both sources, data from VERITAS were analyzed for the VHE regime. 
The study of 3C 66A extends from 2007 to 2015, resulting in a 12 standard deviation (a) detec¬ 
tion from ~61 observing hours. Analysis of W Comae from 2010 to 2014, totaling ~39 hours, 
resulted in a 6(7 low-flux-state detection. We report on the results from these VHE analyses and 
describe contemporaneous multi wavelength data to be used in further analyses. We comment on 
how these low-flux-state IBL detections fit within the context of the blazar paradigm. 
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1. Introduction 

Blazars are a type of active galactic nuclei (AGN) with collimated relativistic jets, specifically 
ones where the jets are aligned with our line of sight. Blazars are some of the most energetic as- 
trophysical objects in the Universe and are found observationally to have a characteristic two-peak 
spectral energy distribution (SED). The lower-frequency peak in photon flux is associated with 
synchrotron radiation resulting from relativistic electrons accelerating in strong magnetic fields. In 
leptonic models, the second peak is thought to be associated with the inverse-Compton (IC) effect, 
with the seed photons being provided by synchrotron radiation (synchrotron self-Compton, SSC), 
external photon fields (external Compton, EC), or some combination of the two e.g. [1]. 

Historically, blazars have been classified as one of two types based on spectral features: flat- 
spectrum radio quasars (ESRQs) or BE Lacertae objects (BE Lacs). ESRQs exhibit strong optical 
emission lines while BE Lacs generally appear to be featureless. ESRQs are generally thought to 
host more efficient accretion disks and stronger jets resulting in lower synchrotron peak frequencies 
(vf^^^), while BE Lacs typically show weaker jets driven by inefficient accretion and cover a range 
of categorized accordingly: high-, intermediate-, or low-frequency-peaked BE Lacs (HBLs, 
IBLs, or LBLs) [2]. However, it has been argued that some of these distinctions are due to sam¬ 
ple selection effects and the wide variety of blazar observables could be accommodated through 
association of BE Lacs and ESRQs with low-excitation (LERGs)/Eanaroff-Riley (ER) I and high- 
excitation (HERGs)/ER II radio galaxies [3, 4]. 

Thus, one of the outstanding problems within blazar observations is the seeming differentia¬ 
tion of blazars into sub-classes as a function of observables such as position of the overall 
source luminosity and the dominance of the IC peak luminosity relative to the synchrotron peak 
("Compton dominance"). These observables can be related to intrinsic physical properties of the 
system such as the maximum particle acceleration energy, the jet power or accretion mode and the 
relative fraction of the external photon fields to the seed photon density. The synchrotron peak is 
observed to fall in a range of frequencies covering nearly five orders of magnitude, from infrared 
to x-ray wavebands and empirically, the Compton dominance is anti-correlated with [2]. In 
the leptonic model paradigm, the ESRQ, LBL, IBL, HBL sequence requires a decreasing external 
photon field component to model the decreasing Compton dominance; ESRQ jets are in the most 
gas enriched environments requiring the highest EC component while HBLs are typically modeled 
satisfactorily with a simple SSC model requiring no EC component [5]. IBLs are transition objects 
and so far have been detected only in flaring states in which SED models have required an EC com¬ 
ponent [6]. If indeed, Compton dominance can be related to underlying physical observables of the 
blazar system, then it is important to understand how the system changes between a quiescent and 
flaring state. 

Thus, to build a more robust understanding of the nature of blazars, it is essential to study 
these transition sources in greater depth. VERITAS has detected many HBLs, but IBL and LBL 
sources are particularly challenging for VHE observatories to detect outside infrequent periods of 
elevated flux. Since the two-peak SED shifts to lower frequencies for IBL/LBLs, the VHE regime 
falls beyond the inverse-Compton peak, leading to very limited observed flux. Therefore, these 
sources require long-term observations to be detected in a low-flux state. Although detections in 
high-flux states (flares) for these sources already exist, it cannot be assumed that IBLs will exhibit 
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the same spectral behavior in a low-flux state; there is already ample evidence of spectral variability 
for HBLs between high- and low-flux states [7, 8]. In particular, it is important to ascertain whether 
the "quiescent" or low-flux-state SED for IBLs requires an external photon field to appropriately 
model the IC peak. This information can be combined with other observables to yield a more com¬ 
prehensive understanding of blazars overall. 

These proceedings describe the results of the VHE analyses of long-term low-flux-state ob¬ 
servations for two IBL sources: 3C 66A and W Comae. Section 2 provides a description of the 
VERITAS observatory and summarizes the past analyses of these two sources in both the high- 
and low-flux states. Section 3 describes highlights of the low-flux-state analysis of 3C 66A, and 
Section 4 details the analysis of W Comae, presenting a low-flux-state spectrum for this source for 
the first time ever. Einally, in Section 5 we offer our conclusions from the analyses thus far and 
describe future work to be done with these sources. 

2. The VHE Detections of 3C 66A and W Comae 

VERITAS is a ground-based array of four imaging atmospheric Cherenkov telescopes (lACTs) 
at the Ered Lawrence Whipple Observatory in southern Arizona. Each telescope is comprised of 
a segmented 12m primary mirror and a camera made of 499 photomultiplier tubes (PMTs). The 
VERITAS telescopes have a 3.5° field of view and are optimized for the energy range 85 GeV - 
30 TeV [9]. VERITAS was constructed in phases; the initial configuration of the four-telescope 
array (2007-2009), called the "old array" (OA), was changed to the "new array" (NA) in the fall of 
2009 when one of the telescopes was relocated to a more optimal position for improved sensitivity. 
Prior to the fall of 2012, the PMTs were replaced, leading to the "upgraded array" (UA), which 
is in current use. The sources discussed in these proceedings were observed at a 0.5° offset from 
the camera center for purposes of simultaneous background estimation [10]. Eor more details on 
VERITAS, see [11]. 



(a) 3C 66A 



188 187 186 185 184 183 

R.A. (J2000) 


(b) W Comae 


Eigure 1: Significance maps (standard cuts) plotted in RA and Dec (degrees) with significance indicated by color bar. 
The white circles indicate regions excluded from the background estimation because they contain stars or gamma-ray 
sources, (a) The analysis for 3C 66A excluded stars brighter than magnitude 7. (b) For W Comae, stars with magnitude 
brighter than 6 were excluded as well as the gamma-ray source lES 1218+304. 
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MJD MJD 

(a) 3C 66A (b) W Comae 

Figure 2: Light curves of low-flux-state analysis for each source, binned by day and plotted against MJD. Observing 
seasons for VERITAS run from September to June. The integral flux (m~^s~^) was calculated above 200 GeV. (a) Data 
for 3C 66A ranges from September 2007 to January 2015. The dotted line represents the average integral flux (>200 
GeV) from the high-flux state (October 2008 flare) of3C 66A. (b) Data for W Comae ranges from January 2010 to May 
2014, dotted line represents the average integral flux (>200 GeV) from the June 2008 flare [6]. 


A VHE detection of 3C 66A (0.3347 < z < 0.41) [12] was first reported by the Crimean As- 
trophysical Observatory in 1998 [13], and a detection of the source by VERITAS was first reported 
in 2008. During an elevated flux state observed for several days at the beginning of October 2008, 
the source reached over 80 " significance, corresponding to an integral flux (>100 GeV) that was 
^10% the flux of the Crab Nebula [14]. VERITAS first reported a low-flux-state detection (5.4<7) 
of 3C 66A in ^17 hours of observations from January 2010 to June 2012 [15]. 

A detection of W Comae (z = 0.102) was first reported in the VHE band by VERITAS after 
an elevated flux state in March 2008. During the two nights of highest activity, W Comae was 
observed at an integral flux (>200 GeV) equivalent to 9% flux from the Crab Nebula [16]. Shortly 
after, in June 2008, W Comae was observed in an even higher flux state (see table 2), with integral 
flux (>200 GeV) corresponding to 25% Crab [6]. VERITAS was also first to report a low-flux-state 
detection of W Comae (5.60") through analysis of ^31 hours of data January 2010 to June 2012 
[15]. In these proceedings, we extend the low-flux-state analyses for both 3C 66A and W Comae. 

3. Highlights of 3C 66A Analysis 

The most recent analysis of 3C 66A, using VERITAS observations from September 2007 to 
January 2015, consists of ^61 hours (livetime) in a low-flux state and ^6 hours (livetime) in a high- 
flux state, after data quality cuts to remove observations with poor weather or hardware issues. The 
high-flux state (October 2008 flare) was distinguished from the low-flux state by showing that the 
probability of randomly seeing an integral flux level significantly higher than the mean flux three 
days in a row is negligible. Analysis of this source utilized a combination of soft (<200 GeV) and 
standard (>200 GeV) cuts, determined by studying data from the Crab Nebula scaled to 10% flux 
as a toy model. Soft cuts allow for events of lower minimum size (as measured by digital counts) 
than the standard cuts and also use a larger angular region around the source to account for poorer 
angular resolution at this lowered threshold. 
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3C 66A Expos. Signif. Normalization C Photon Index Curvature Integral flux 

(hrs) (7 (m“^s“^TeV“^) F /3 

(a) Low 61 12 (3.62±0.33) X 10-^ 4.67±0.24 3.65±1.24 (7.08±0.74) x 10"^ 

(b) High 6 13 (1.40±0.16) X 10-^ 5.18±0.34 2.29±1.72 (2.98±0.37) x lO"’ 


Table 1: Statistics of low- and high-flux states of3C 66A. Exposure is the total livetime in hours after data quality cuts. 
Analysis for both states uses a combination of soft cuts (<200 GeV) and standard cuts (>200 GeV). A log-parabola 
function was used to model the data. The lower limit for the integral flux is 200 GeV; the integral flux corresponds to 
^3% Crab for the low-flux state and r^l2% Crab for the high-flux state. All errors shown are statistical. 


A position-based analysis of VERITAS observations has shown that the center of the excess is 
associated with 3C 66A, ruling out radio galaxy 3C 66B at 95% confidence. Analysis of the low- 
flux state resulted in a \la detection. The photon spectrum for the low-flux state was found to be 
best fit with a log-parabola function of the form dN/dE = C{E where C is the 
normalization constant, Eq is the decorrelation energy, F is the photon index, and j8 is the curvature 
parameter. The decorrelation energy, Eq = 0.27 TeV in this analysis, represents the energy at 
which the calculated flux is least affected by the choice of spectral index. See Table 1 for the other 
parameters and Figure 3a for the low-flux-state spectrum. Since a log-parabola model is preferred 
over a power law at 4.1(7, we conclude that there is significant curvature to the photon spectrum. 
From the results of the fits, we also note that the photon index F and curvature parameter j8 for 
the low- and high-flux states agree within statistical uncertainty, showing that there is currently no 
observable change of shape in the VHE spectrum between states. A more detailed description of 
the analysis of 3C 66A and conclusions drawn will appear in an upcoming publication from the 
VERITAS collaboration. 

4. Detailed VHE Analysis of W Comae in Low-Flux State 

Beyond the previous low-flux-state detection of W Comae in 2012, this analysis includes two 
seasons of observations (2012/13 and 2013/14) from the upgraded VERITAS array, resulting in 
^39 hours of total livetime after data quality cuts. All observations by VERITAS of W Comae 
over the span of 2010-2014 are defined as the low-flux state, which was determined by comparing 
the integral flux of the high-flux state (June 2008 flare) to the daily light curve from 2010 to 2014 
(see Figure 2a). 

Two independent analyses of W Comae were performed using standard cuts and produced 
consistent results. The VHE significance map (Figure lb) shows the VERITAS field a view in 
the vicinity of W Comae. A reflected region background model was used, where background is 
estimated from circular regions of equal size and at equal distance from the camera center as the 
source (0.5° offset) [10]. From this analysis, observations yielded 164 excess events (Non = 883 
and Noff = 6933, average normalization factor a = 0.1014), resulting in a 6.2(7 detection and in¬ 
tegral flux (>200 GeV) of (2.81 ±0.53) x 10“^ m“^s“^ from ^39 hours of livetime. The first ever 
low-flux-state VHE spectrum of W Comae is presented in Figure 3b. 

The spectrum was fit with a power law model of the form dN/dE — C(E/Ef)~^, where C is 
the normalization constant, Eq is the decorrelation energy (Eq = 0.32 TeV in this case), and F is 
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Figure 3: Low-flux-state VHE spectra ofW Comae and 3C 66A, shown in the E^dN/dE representation (TeVm~^s~^). 
Bins with significance >2a and at least 5 signal events after background subtraction are plotted. All errors plotted are 
statistical, (a) 3C 66A is detected in a low-flux state at 12g using soft (<200 GeV) and standard (>200 GeV) cuts with 
61 hours of observations. A log-parabola function is used to model the data, (b) W Comae is detected in low-flux state 
at 6.2 G using standard cuts with 39 hours of livetime. A power law function is used to model the data. 


the photon index. The best fit of the data using this model results in C = (6.29 ± 1.41^(a() x 10 ^ 
m“^s“*TeV“* andT = 3.\±0.6stat- We note that the value of the photon index in the low-flux state 
for W Comae is in agreement with the photon index for the June 2008 flare (F = 3.68 ±0.22^^^^) [6] 
within statistical uncertainty, showing that we currently find no significant change in the spectral 
shape between states in the VHE waveband. 

Future work with W Comae will continue to provide greater insight into the underlying physics 
behind IBL sources. This includes extending the analysis both within the VHE waveband and to 
other wavelengths. In the VHE waveband, we intend to integrate data from the old array (2007/08 
and 2008/09 seasons) and also the most recent season of data from the upgraded array (2014/15). 
Analysis of the most recent season is already underway; our preliminary results show W Comae to 
be unusually faint in VHE during 2014/15. 

In addition to continued work with VERITAS data, we are in the process of gathering publicly 
available multiwavelength data from radio to HE gamma rays for the purpose of constructing the 


W Comae 

Exposure 

Significance 

Normalization C 

Photon Index 

Integral Flux 


(hrs) 

G 

(m-^s-'TeV-') 

r 


(a) Low 

39 

6 

(6.29 ±1.41) X 10“^ 

3.1 ±0.6 

(2.81 ±0.53) X 10“^ 

(b) High 

4 

10 

(6.5 ±0.9) X 10“’ 

3.68 ±0.22 

(5.7 ±0.6) X 10-’ 


Table 2: Statistics of low-flux and high-flux (June 2008 flare) states ofW Comae. Eq = 0.4 TeV for June 2008 flare 
[6]. The data is fit with a power law model for both states. The lower limit for integral flux is 200 GeV; the integral flux 
corresponds to Crab for the low-flux state and r^25% Crab for the high-flux state. All errors shown are statistical. 
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full multi-band SED most representative of the low-flux state. In the HE regime, we have already 
begun an analysis of W Comae using data from the Fermi Large Area Telescope (EAT). As a Fermi 
source of interest, W Comae is one of many blazars also monitored by Swift in the x-ray through 
optical wavebands {Swift XRT and UVOT). Upon initial study. Swift public x-ray data shows a 
decline in integral flux over the course of the most recent season of VERITAS data (January-April 
2015), matching our preliminary VHE results. Several optical observatories have long-term data on 
W Comae, including the Catalina Sky Survey, the ELWO 48” Telescope, and Tuorla, and a combi¬ 
nation of these datasets will provide contemporaneous observations in the optical regime. In radio, 
W Comae is observed by the Owens Valley Radio Observatory (OVRO) at 15 GHz approximately 
twice each week, and thus we expect there to be contemporaneous data with the VHE observa¬ 
tions. We will continue to investigate this dataset and query other radio observatories for more data 
to include in a low-flux-state SED. Study of this full SED will provide insight into the spectral 
variability of W Comae between high- and low-flux states and also can lead to understanding of 
the intrinsic jet strength. 


5. Conclusions 

Through long-term, low-flux-state analyses of 3C 66A and W Comae, we contribute the VHE 
portion of a future SED for each of these IBL sources in their low-flux states. VERITAS obser¬ 
vations of 3C 66A totaling ^61 hours of live time in the low-flux state led to a 12<7 detection. A 
log-parabola model was found to better fit the spectrum than a power law model, showing signifi¬ 
cant curvature in the VHE spectrum. Observations of W Comae totaling ^39 hours led to a 6.2(7 
detection, allowing for a VHE spectrum in the low-flux state to be derived for the first time. We see 
no significant curvature in the VHE spectrum for W Comae, with the data being best fit to a power 
law model. Eor both sources, there is no significant difference in the spectral shape between the 
low- and high-flux states in VHE. Continued observations of these sources in a low-flux state will 
allow for greater resolution in the VHE spectra and better constrained parameters. 

These results provide the initial step towards multiwavelength analyses. The VHE spectrum 
in combination with other wavelengths will define the shape and locations of both the synchrotron 
and inverse-Compton peaks, providing insight on the nature of IBLs in multiple ways. It has been 
shown that there is a shift in the SED of HBL sources during flaring states, and so we can deter¬ 
mine whether this effect also exists for IBLs once a low-flux-state SED is produced. Eurthermore, 
the relative luminosities of the synchrotron and inverse-Compton peaks can offer information on 
the strength of external photon fields, which could provide evidence for whether or not IBLs are 
driven by weak or strong jets. The study of IBLs as transition objects in synchrotron peak frequency 
provides an essential step towards developing greater understanding of the overall blazar paradigm. 
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